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Summary 

Data partitioning and modified stepwise regres- 
sion were applied to recorded flight data from a Royal 
Aerospace Establishment (RAE) high-incidence re- 
search model (HIRM). An aerodynamic model struc- 
ture and the corresponding stability and control 
derivatives were determined for angles of attack 
between 18° and 30°. Several nonlinearities in 
angles of attack and sideslip as well as a unique roll- 
dominated set of lateral modes were found. All flight- 
estimated values were compared with available wind- 
tunnel measurements. 

Introduction 

The Royal Aerospace Establishment (RAE) is 
currently leading an extensive program of research 
into the high-angle-of-attack aerodynamic behavior 
of current and future combat aircraft configurations 
(ref. 1). The main thrust of the program is to pro- 
duce an adequate mathematical model of the aircraft 
aerodynamics at high angles of attack so that control 
laws can be designed for enhanced maneuverability 
and departure prevention. The program includes the- 
oretical studies (ref. 2), extensive wind-tunnel test- 
ing (refs. 3-6), and flight testing of unpowered drop 
models (ref. 7). One of the vehicles for investiga- 
tion is the high-incidence research model (HIRM), a 
canard/swept back- wing configuration. 

While theoretical analysis, wind-tunnel testing, 
and short test flights (20 sec or less) are being carried 
out in the United Kingdom (UK) by the RAE, the 
responsibility for extended test flights (up to 160 sec) 
and flight data analysis is shared with the National 
Aeronautics and Space Administration (NASA) at its 
Dryden Flight Research Facility (DFRF) and Lang- 
ley Research Center (LaRC). Flight tests were car- 
ried out in 1983 and 1986 under Task 18 of the 
NASA/RAE Joint Aeronautical Program. As a re- 
sult of unpredicted, violent wing rock episodes in 
several of the 1983 flights, Task 32 was added to 
the joint program for the purpose of determining the 
cause of that wing rock using LaRC’s system iden- 
tification expertise. DFRF provided flight test sup- 
port, and LaRC planned the 1986 flight maneuvers 
and applied its system identification expertise to the 
acquired flight data from 1983 and 1986 for mathe- 
matical modeling. Some results from the analysis of 
the 1983 series were published previously (ref. 2) and 
were used to plan the 1986 test program. 

The purpose of this report is to present results 
from the analysis of flight data from both the 1983 
and 1986 series. The aim of that analysis was to 
determine an adequate model of the aerodynamics 
of the HIRM vehicle in the tested flight regimes. 


This model should explain all analyzed data and it 
should also have good prediction capability for data 
not included in the analysis. In order to obtain an 
adequate model, three procedural steps were taken: 
(1) data were partitioned into subsets, each com- 
prising a small angle-of-attack range; (2) stepwise 
regression was applied to each subset to determine 
model structure and stability and control derivatives; 
and (3) additional research flights were designed and 
flown, and steps 1 and 2 were applied to them. The 
results are mainly in the form of body-axis stability 
and control derivatives. This report will include a 
sample of flight data time histories, a discussion of 
the data analysis techniques, and all the results ob- 
tained from applying those techniques to the flight 
data. All results will be compared with available 
wind-tunnel data. The report is organized as fol- 
lows: After this introduction, the flight vehicle, flight 
tests, and flight data will be discussed. Next, avail- 
able wind-tunnel data will be summarized and flight 
data analysis methods will be presented. The main 
results of the data analysis follow, and a concluding 
section summarizes the results and provides sugges- 
tions for further work. 

Symbols 

Values are given in SI Units, but they are occa- 
sionally given in U.S. Customary Units where con- 
sidered useful. Measurements and calculations were 
made in SI Units. 


&x i dyi 

longitudinal, lateral, and verti- 
cal accelerations, respectively, 
g units 

b 

wing span, m 

C a 

general aerodynamic force and 
moment coefficient 

Cl 

lift coefficient, L/qS 

Ci 

rolling-moment coefficient, 
M x /qSb 

Cm 

pitching- moment coefficient, 
My/qSc 

C n 

yawing- moment coefficient, 
Mz/qSb 

Cx 

longitudinal- force coefficient, 

Fx/qS 

Cy 

lateral-force coefficient, Fy/qS 

Cz 

vertical- force coefficient, 

F Z /qS 

c 

wing mean aerodynamic chord, 
m 



F 

statistical F-values 

f x ,f y ,f z 

forces along longitudinal, 
lateral, and vertical body axes, 
respectively, N 

/( ) 

function 

9 

acceleration due to gravity 
(Ig % 9.81 m/sec 2 ), m/sec 2 

txJyJz 

moments of inertia about lon- 
gitudinal, lateral, and vertical 
body axes, respectively, kg-m 2 

Ixz 

product of inertia, kg-m 2 

i 

= v^T 

J 

cost function 


damping factors for Dutch 
roll, roll, and spiral modes, 
respectively 

k 1 .k 2 .k 3 

damping factor for three 
lateral modes found in flight 
data 

L 

lift force, N 

M X , My , 

rolling, pitching, and yawing 
moments, N-m 

m 

mass, kg 

N 

number of data points 

N Re 

Reynolds number based on 
wing reference chord 

P 

roll rate, rad/sec or deg/sec 

9 

pitch rate, rad/sec or deg/sec 

9 

dynamic pressure, \pV 2 , Pa 

R 2 

square of multiple-correlation 
coefficient 

r 

yaw rate, rad/ sec or deg/sec 

S 

wing area, m 2 

s 2 

residual mean square 

t 

time, sec 

U, V, w 

longitudinal, lateral, and 
vertical airspeed components, 
respectively, m/sec 

V 

airspeed, m/sec 

X 

matrix of independent 
variables 


x,v,z 

longitudinal, lateral, and 
vertical body axes, respectively 

Y 

vector of measurements of 
dependent variable 

y 

dependent variable in regres- 
sion equation 

a 

angle of attack, rad or deg 

0 

sideslip angle, rad or deg 

Sc 

symmetric canard deflection, 
rad or deg 

^control 

any control variable, rad or 
deg 

^cd 

differential canard deflection, 
rad or deg 

Sd 

differential tail deflection, rad 
or deg 

Sh 

symmetric tail deflection, rad 
or deg 

Sr 

rudder deflection, rad or deg 

e 

equation error vector 

0 

unknown parameter vector 

Oi 

stability and control derivative 

6 

pitch angle, rad 

A 

eigenvalue 

P 

air density, kg/m 3 

V 

residual (difference between 
values of measured and pre- 
dicted or computed variable) 


roll angle, rad 


heading angle, rad 

U) 

frequency of lateral oscillatory 
mode from flight data, 1/sec 

U)£> 

frequency of Dutch roll mode, 
1/sec 

Abbreviations: 

BAe 

British Aerospace 

c.g. 

center of gravity 

HIRM 

high-incidence research model 

MSR 

modified stepwise regression 

RAE 

Royal Aerospace Establish- 
ment (formerly Royal Aircraft 
Establishment) 
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SAS stability augmentation system 

SR stepwise regression 

A dot (•) over a symbol denotes a derivative with 
respect to time. A bar (') over a symbol indicates an 
average value. A circumflex (*) denotes an estimated 
value. 

Derivatives of the aerodynamic coefficient C a 
(where a = X, Y, Z, m, /, and n) referenced to 
a system of body axes with the origin at the airplane 
center of gravity are given as follows: 



dC a 


_ 8C a 

Ca p — 


Ca q : 
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Ca a 
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Flight Vehicle and Flight Test 

HIRM is a three-surface, close-coupled canard/ 
sweptback-wing drop model representative of a class 
of future fighter aircraft configurations. The canard/ 
swept-wing combination allows for superior transonic 
performance, and an all-moving tail provides en- 
hanced pitch control. Drawings of the vehicle are 
presented in figure 1. Geometric, mass, and inertia 
characteristics (where the c.g. is at 12.5-percent c) 
are given in the following table: 


Wing area, m 2 2.062 

Wing span, m 2.504 

Wing mean aerodynamic chord, m .... 0.868 

Mass, kg 248.8 

lx, kg-m 2 17.8 

ly, kg-m 2 151 

l z , kg-m 2 164 

Ixz > kg-m 2 -0.19 


Longitudinal control is achieved by a symmetric 
canard and all-moving tail (used either separately 
or together). Lateral-directional control is achieved 
through the use of a differential canard, differential 
horizontal tail, and rudder. The flight vehicle was 
sized so that a full-scale model could be tested in the 
wind tunnel at flight Reynolds numbers. 

Flight instrumentation included three rate gyros, 
roll and pitch attitude gyros, three linear accelerom- 
eters, angle-of-attack and sideslip sensors, and a 
pitot static probe; potentiometers were installed at 


each control surface to measure surface deflections. 
Flight data were telemetered to a ground station 
and recorded for future processing. There was no 
ground- to- vehicle radio uplink. Therefore, all flights 
had to be carefully planned using the best available 
mathematical model to design a preprogrammed set 
of control inputs to achieve desired flight test condi- 
tions and responses. After a flight had been planned, 
trim conditions and simple responses were checked 
in short flights (approximately 20-sec duration) at 
RAE Larkhill in the UK. The major test flights took 
place in two series (1983 and 1986) at the NASA Dry- 
den Flight Research Facility (DFRF) in the United 
States. The flights at DFRF normally lasted between 
120 and 140 sec each. 

The 1983 series investigated the response of 
HIRM to differential canard and differential tail dou- 
blets and rudder pulses for the two mean canard set- 
tings of 0° and —10° and to mean tail settings to 
trim the vehicle at angles of attack greater than 20°. 
During several of these flights, a large wing rock type 
of lateral oscillation occurred spontaneously. 

The 1986 series was planned to investigate this 
wing rock region further as well as test several control 
laws for departure prevention and for angle-of-attack 
and bank-angle control. The control laws for depar- 
ture prevention (DEPS) were formulated to prevent 
departure from controlled flight by limiting angle of 
attack. The other set of control laws (HIRM aero- 
dynamic parameter identification, or HAPI) was de- 
signed to maintain steady longitudinal flight during 
lateral responses to yaw and roll control inputs and 
to prevent excessive bank angles (greater than 30°). 
Canard settings for the 1986 series were mostly 0° 
and —5° with some data taken for 6 C = —3°. There 
were no differential canard inputs designed for system 
identification; the rudder and differential tail were 
employed only to control excessive bank angle. Sym- 
metric tail inputs were limited to steps required to 
change the trim angle of attack. 

A typical flight at DFRF began by towing HIRM 
by helicopter to an altitude of approximately 
10 000 ft. The vehicle was then released to glide 
back to Earth while executing its preprogrammed 
set of control instructions along with those inputs 
demanded by the augmentation system. At approx- 
imately 1000 ft above ground level, a parachute was 
deployed to slow the vehicle and lower it to the 
ground. Air bags were deployed beneath the vehi- 
cle to cushion its landing. 

Flight Data 

Recorded flight data were sent from the 
Dryden Flight Research Facility to the RAE (at 
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Farnborough, UK) where they were converted to en- 
gineering units and digitized to a time between sam- 
ples of 0.0116 sec. Those digitized data were then 
sent to NASA LaRC on several 9-track magnetic 
tapes for data analysis. Time histories for selected in- 
put and response variables are given in figure 2 for the 
first 35 sec of one of the 1983 flights. The rapid pitch- 
up and settling out to the trim angle of attack after 
release from tow (at approximately 2.6 sec) is seen. 
At approximately 10 sec into the flight, an exagger- 
ated Dutch roll or wing rock type of motion starts 
to build. The oscillatory motion is disturbed only 
briefly by the sharp differential tail doublet (origi- 
nally planned to excite rolling motion for parameter 
identification) at 23 sec into the flight. For com- 
parison, time histories from the first 40 sec of one 
of the augmented 1986 flights are given in figure 3. 
Again, there is a large longitudinal response to the 
release from the tow line. However, the active aug- 
mentation is visible as increased control activity — 
particularly in differential tail and rudder — to keep 
the bank angle controlled when changing trim con- 
ditions. The continuation of the 1983 flight is given 
in figure 4 (with a new scale) in which a wing rock 
spontaneously develops as the trim angle of attack 
reaches down into the 24° region. The end time of 
figure 3 and the start time of figure 4 appear to over- 
lap but do not, because figure 3 was shifted 2.6 sec 
for the purpose of plotting. 

The data used in the analysis presented in this 
paper are drawn from five 1983 flights (RAE des- 
ignations HD1, HD3, HD5, HD2, and HD4) and 
from four 1986 flights (HD8, HD 19, HAD9, and 
HAD16). Flights HD1, HD3, and HD5 have -10° 
canard settings, and flights HD2 and HD4 have a 0° 
canard setting. For the 1986 series all flights included 
both 0° and -5° canard settings. Flights HD8 and 
HD19 were unaugmented, whereas flights HAD9 and 
HAD 16 were augmented. 

The 1983 trials contained several well-planned 
system identification maneuvers in which the vehicle 
was excited first by the release from the helicopter 
and then from trim using programmed control in- 
puts. Those control inputs were elevator square-wave 
doublets designed to excite the short-period longitu- 
dinal dynamics of the vehicle and rudder/differential 
tail square-wave doublet combinations to excite lat- 
eral motion. An individual maneuver comprised the 
input doublet(s) and response to that input. Each of 
the individual maneuvers was analyzed by itself. The 
1986 trials, designed to gently and carefully probe 
the wing rock regions, had only preprogrammed step 
inputs of less than 1° from the horizontal stabilator 
(longitudinal), whereas lateral inputs were initiated 
only by the flight control system computer as the 


bank angle approached ±30° or the trim angle of 
attack was changed. An example of an individual 
maneuver from the 1983 series is the differential hor- 
izontal stabilator doublet shown in figure 5(a). Fig- 
ure 5(b) shows the type of control surface motion 
that had to be utilized as initiating an individual 
maneuver in the 1986 test data. 

Wind-Tunnel Tests 

The RAE conducted extensive wind-tunnel test- 
ing of the HIRM. Data were taken for 0.22-, 
0.36-, 0.44-, and full-scale models. Early static tests 
were conducted on a 0.36-scale preliminary version 
of HIRM in the RAE (Farnborough) No. 2 ll^- by 
8-Foot Tunnel (0° < a < 40°) and in the RAE 
(Bedford) 13- by 9-Foot Tunnel (0° < a < 90°). The 
purpose of these tests was to define a final HIRM ge- 
ometry and to develop baseline data on lift, drag, and 
control effectiveness especially with regard to trim 
angles of attack. These baseline data were then used 
to define the test matrix for the definitive HIRM ge- 
ometry. Data were taken at 2° increments in angle 
of attack for 0° < a < 40° and at 5° increments for 
40° < a < 90°. Based on these tests, a final geome- 
try for HIRM was defined and a 0.44-scale model was 
built. Before the full-scale free flight tests, this model 
underwent static and oscillatory tests in the RAE 
(Bedford) 13- by 9-Foot Tunnel and rotary tests 
in the BAe (Warton) 18-Foot Tunnel. After the 
first series of free flights at DFRF, a 0.22-scale 
model was tested on the new rotary rig at the RAE 
(Bedford) 13- by 9-Foot Tunnel. Tunnel wind speed 
was 70 m/sec allowing for N^ e = 1.5 x 10 6 (based 
on c) which is close to the full-scale model flying at 
40 m/sec at an altitude between 14 000 and 2000 ft 
(with Nft e = 1.7 x 10 6 to 2.3 x 10 6 ). Additional 
rolling- and yawing- moment coefficients were mea- 
sured for the full-scale free flight model in the RAE 
(Farnborough) 16.5-Foot Tunnel. A summary of 
available derivatives, their sources, and the figure in 
which each is used is given in table I. 

Analysis Methods 

A general flow chart for system identification 
analysis is given in figure 6. 

First, an experiment is designed using all available 
a priori knowledge of the physical system to be 
identified. This knowledge includes results from 
previous experiments. Once the experiment has been 
performed and data have been gathered, the data are 
run through a compatibility check. The purpose of 
the compatibility check is to determine whether all 
data channels are self and cross consistent. The check 
should also allow for the estimation of bias or scale 
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factor errors in the data so as the make the measured 
data correctable. 

After the data are brought into the most con- 
sistent shape possible, the analysis continues with 
model structure determination and parameter esti- 
mation. Here, a priori knowledge of the physical 
system and model building software such as step- 
wise regression are used to determine a set of models 
that are qualitatively adequate to explain the behav- 
ior of the system as recorded in the measured data. 
The model is completed by estimating any constants 
(parameters) associated with the determined model 
structure. This model is then verified by assessing 
its prediction capability using a data set that was 
not employed in the identification. If the verification 
is good, the process is complete, with the model be- 
ing that which was verified. If the verification result 
is lacking, then the analyst must backtrack either to 
the experiment design stage or simply to the model 
structure determination stage. The return to the lat- 
ter stage can pass through a data restructuring stage 
(such as data partitioning) in which the measured 
data are separated or combined in innovative ways to 
allow for better model structure determination. This 
inner loop continues until a successful verification is 
achieved or until all experiments and data restruc- 
turing have been exhausted. 

The following paragraphs will focus on a data 
compatibility check; a model structure determination 
and parameter estimation; and, as was necessary for 
the analysis of these data, partitioning. 

Data Compatibility 

Each flight from the 1983 and 1986 series was first 
subjected to a test for data self-consistency (data 
compatibility). This check is performed by integrat- 
ing the aircraft kinematic equations using measured 
values of angular rates and linear accelerations as in- 
puts that are integrated to compute output variables 
consisting of Euler angles, total airspeed, angle of 
attack, and angle of sideslip. An algorithm reported 
in reference 8 was used. In addition to the simple 
integration, the algorithm allows the option of esti- 
mating values for bias and scale-factor errors on the 
input and output data that would minimize the sum 
of squares of differences between measured and com- 
puted values (residuals). Along with identifying bias 
and scale-factor corrections, the data compatibility 
check allows for an assessment of gross data prob- 
lems such as sign reversals. If several sections of a 
flight yield consistent estimates of bias or scale fac- 
tors, the data from that flight are corrected by ap- 
plying that bias or scale-factor correction to the data 
prior to analysis. After these data were corrected, the 
analysis proceeded with the determination of model 


structure and the estimation of corresponding stabil- 
ity and control derivatives. 

Data Restructuring 

The general form of the aircraft equations of mo- 
tion and aerodynamic model equations is found in 
the appendix. Nonlinear dependence of the aerody- 
namic force or moment coefficients on one or more 
variables leads to difficulty in postulation of a proper 
model structure, whereas lack of good information 
or lack of harmonic content in a maneuver leads to 
identifiability problems with some parameters. The 
large amplitude motions that often characterize drop 
model flight data can pose both of these problems. 
Multiple nonlinear dependence causes the entry of 
nonlinear terms in more than one variable when the 
modified stepwise regression (explained below) is ap- 
plied to data from a maneuver. The variation of any 
potential independent variable (such as a, (3, etc.) in 
a maneuver can be reduced for the purpose of anal- 
ysis by dividing the maneuver into several subsets 
(with each subset covering a fraction of the total ma- 
neuver range) and then analyzing each subset inde- 
pendently. This process is called partitioning. 

To understand the basis for partitioning, let 


y(0 = /[*i(0, *2 (t), ..., *„(<)] (1) 


where x* to x n are the regressors formed from the air- 
plane response and control variables. Now suppose 
that for any time t we want to eliminate the depen- 
dence of y(t) on one of the Xj(t), say x p (t). Then, 
by partitioning we mean redefining y(t) on proper 
subsets of {xi, X 2 , x n } as 


y(x p ,t) 


r h [*i(0i *2(0, *p-i(0, * P + i(0. •••».*«(<)] 

(*po < x p < x Pl ) 
h [*l(0. *2(0 *p-l(0< *p+l(0 *«(0] 

(jCpi < Xp <■ XP2) 


fm [*l(*h x iW x p-l (*)' *p+l(0i •••• *n(0] 

, ( x Pm-l < X P< x Pm) 

( 2 ) 

where 


*> = =*±1 p* (i = 1,2, ...,m) 

That is, each (n + l)-tuple in (xj,j/) is reduced 
to several n-tuples — each associated with a particu- 
lar value on a small range of x p . The supposition 
is that as the range defined by x Pi+l - x Pi becomes 
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smaller, the variation in / due to x p becomes less sig- 
nificant and reaches a level where it can be neglected* 
For example, an aircraft might perform a mostly lat- 
eral maneuver but with angle of attack a varying 
between 20° and 30°. Because of separation effects 
on the lifting surfaces in this a region, we could ex- 
pect that the lateral aerodynamic force and moment 
coefficients Cy,C\, and C n might well depend on a 
in a nonlinear way, i.e., 


Cn Cn(a, /3, p, r, ^control) (3) 

Then, to partition, one would simply analyze the 
data in separate groupings as follows: 


(CVi)f*=21° — Cn(A Pi ^control) 
(Cri)a=23° — Cri(3< p, r, <5 con trol) 


(20° < a < 22°) 'j 


(22° < a < 24°) I 


(4) 


(C7*)a=29° — P * r . ^control) (28° < a < 30°) ) 


After partitioning for some individual maneuvers, 
the information content of the data to be analyzed 
may be inadequate. Indicators of this problem are a 
lack of identifiability for some parameters (because 
of large scatter or large standard error of parameter 
estimates) or cross plots of variables that indicate a 
lack of excitation in one or more variables. An ex- 
ample of cross plots of roll rate and sideslip angle 
(p versus (3) for eight bins after partitioning one of 
the 1986 flights is given in figure 8. Figure 8(a) shows 
data that are minimally adequate in distribution and 
range for identification. Figures 8(d) and 8(e) are en- 
tirely adequate in both range and distribution within 
that range. Figure 8(h) is adequate in distribution 
pattern but is limited in range for both f3 and p. To 
increase the information content of the data being 
analyzed, several maneuvers can be combined and 
the combined set be partitioned; or, equivalently, the 
partitioned sets from equivalent a bins from several 
flights can be combined for model structure deter- 
mination and parameter estimation using modified 
stepwise regression. 


That is, all data corresponding to 20° < a < 22° are 
put into one group for analysis, data corresponding 
to parts of the maneuver in which 22° < a < 24° are 
put into a second group, and so forth, until all data 
have been accounted for. 

If any grouping still appears to be dependent on 
a, it can be subdivided further (assuming that a suf- 
ficient number of data points exist). In this exam- 
ple one can now analyze C n = f(/3, p, r, 6 con trol) at 
characteristic values of a given by the mean value of 
a for each grouping. To be more specific, consider 
an application of the above process to data with the 
angle-of-attack time history shown in figure 7. All 
data corresponding to angles of attack less than 20° 
are put in bin 1, data with 20° < a < 22° are put in 
bin 2, and so forth, until all data with a > 30° are 
put in bin 7. Notice that this partitioning procedure 
will often lead to non-time-contiguous data being an- 
alyzed in the same bin (as with bin 2 for sections of 
data at a time approximately equal to 3, 5, 19, 22, 
and 25 sec). 

After the bins or subsets are established, the 
model structure determination and parameter esti- 
mation can proceed. The model structure determi- 
nation is still necessary since large variations also 
may have occurred in variables other than the one 
on which the partitioning is based (for example, large 
variations in 0 while partitioning is with respect to 
a). The model structure determination and param- 
eter estimation then proceed by applying a modified 
stepwise regression (MSR) algorithm to each bin or 
subset of partitioned data (ref. 9). 


Model Structure Determination and Parameter 
Estimation 

All the data from each bin are analyzed using a 
modified stepwise regression (MSR). As a modified 
version of the linear regression, this method can de- 
termine the structure of an aerodynamic model and 
estimate the model parameters. The determination 
of an adequate model (a model that fits the data and 
has good prediction capabilities) for the aerodynamic 
coefficients includes three steps: the postulation of 
terms that might enter the model, the selection of 
a model, and the verification of the model selected. 
The general form of aerodynamic model equations 
can be written as 

y(t) = 0 o + 0i xi (t) ... +9 n x n (t) (5) 

where y(t) represents the resultant coefficient of aero- 
dynamic force or moment. In the polynomial rep- 
resentation of the aerodynamic coefficient, 6\ to 
On are the stability and control derivatives. (The 
linear stability and control derivatives are of the 
form C Axi « dCj\/dxi where A — X,Y, 
or n and — a, /?, p, r, or <5 con t ro p nonlinear 
derivatives are of the form d ni C A /dx\ l ...dx^ 1 where 
2, 3 , ... , riy i, ..., k — 1 , 2, ..., n; ij, ..., h\ — 

1,2, 3, 4, 5; and n\ — i\ + ... + k\. For example, 
in d 2 Ci/dad(3,C A = Q,ni = 2;x\ = a,x 2 - f); 
and ii = 1,* 2 = L The symbol Oq is the value of 
any particular coefficient corresponding to the refer- 
ence flight conditions, and X{ to x n are the regressors 
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formed from the airplane output and control vari- 
ables and their combinations. 

Postulating the aerodynamic model equations by 
selecting the candidate regressors is followed by 
the determination of significant terms among the 
candidate variables and the estimation of the corre- 
sponding parameters. The variable chosen for entry 
into the regression equation is the one that has the 
largest correlation with y after adjusting for the effect 
on y of the variables already selected. The parame- 
ters are estimated by minimizing the cost function 


N 


J SR = /I 
i = 1 



( 6 ) 


where N is the number of data points and t + 1 is the 
number of parameters in the regression equation. 

At every step of the regression, the variables in- 
corporated into the model in previous stages and 
a new variable entering the model are reexamined. 
Any variable that provides a nonsignificant contri- 
bution as given by its F-statistic (due to correlation 
with more recently added terms) is removed from the 
model by the algorithm. The process of selecting and 
checking variables continues until no more variables 
are admitted to the equation and no more are re- 
jected. Experience shows, however, that the model 
based only on the significance of individual param- 
eters in the model in equation (6) can still include 
too many terms and therefore may have poor pre- 
diction capabilities. Therefore, three quantities are 
examined for each model structure as criteria for an 
adequate model, and they are presented as follows: 


1. The computed value of the F-statistic. This is 
given as the ratio of regression mean square to resid- 
ual mean square and is calculated as 


F = 



(7) 


where 


y = jjlt v(*) 

i = 1 

* 2 = a rb; £ [m - vw 2 

i=l 


and the superscript T denotes a transpose matrix. 
Equation (7) was developed for the model given by 


equation (1) and can be expressed in matrix notation 
as 


Y = XG + e (8) 


The F-value can be associated with an estimate of 
the standardized total mean-square error of estima- 
tion for the data X as shown in reference 10. The 
model with a maximum F-value has already been 
recommended in reference 11 as the “best” one for a 
given set of data. 


2. The value of the square of the multiple- 
correlation coefficient (R 2 ). This indicates the pro- 
portion of the variation due to terms other than Oq 
in the model. The value of R 2 is obtained from the 
equation 

0 t X t Y - Ny 2 


R 2 = 


YT’Y - Nil 2 


(9) 


3. The residuals e(i). For an adequate model, 
the sequence of the residuals e(i) should be close to a 
random sequence that is uncorrelated and Gaussian. 


The use of these three criteria is demonstrated in 
references 12 and 13. 


Results 


Data Compatibility 

Data received from the United Kingdom exhib- 
ited good compatibility for both the 1983 and 1986 
flight series. Because of the large excitation just after 
release from the helicopter, due to the mismatch of 
tow airspeed and angle of attack with those for free 
flight, the beginning of each flight offered the best 
opportunity for a compatibility check. An example 
of the highly excited beginning of one of the 1986 
flights is given in figure 9. Figure 9(a) gives the mea- 
sured and computed time histories for output vari- 
ables V, a, 9, (3, and 4>. Figure 9(b) gives the residual 
time histories (the difference between the measured 
and computed values for the corresponding output 
variable in fig. 9(a)) for V,a,6,(3, and 4>. Figure 9(c) 
gives the time histories of the linear accelerations and 
angular rates that are the inputs to the data compati- 
bility algorithm. The computed values in figure 9 are 
for the case in which biases were estimated and re- 
moved from a x ,a y ,a z ,p,q, and r. The estimated bias 
values for a x ,a.y,a z ,p,q, and r for this part of the 
flight were O.OI83, — 0.015<7, — 0.089<7, 0.02 rad/sec, 
-0.002 rad/sec, and 0.03 rad/sec, respectively. 
These bias values are small relative to the ranges of 
the analyzed maneuver. Moreover, when several dif- 
ferent flights or different sections of data from one 
flight were analyzed, no consistent bias values could 
be found. Therefore, it was determined that no bias 


7 



corrections would be made before analysis. The ef- 
fect of this decision on the accuracy of estimated sta- 
bility and control derivatives is thought to be small 
(heuristically) but cannot actually be determined. 

Aerodynamic Model Structure 

Lateral . Model structures were determined by 
applying MSR to sections of data identified as “in- 
dividual maneuvers 1 ’ and to bins of data created by 
partitioning. The individual maneuvers were chosen 
to be portions of a flight in which there were sig- 
nificant control surface movements (±2° or 3°) and 
significant vehicle response to those movements (sig- 
nal to noise greater than 1 visually on plotted time 
histories). An adequate width of bins for the parti- 
tioned data was determined, in part, by the results of 
several MSR applications described as follows: First, 
a nominal bin width of 2° angle of attack was chosen. 
This 2° width often precipitated the inclusion of non- 
linear terms that indicated angle-of-attack depen- 
dence (such as pa, pa 2 , 0a, and so forth). Then, the 
bin width was reduced and MSR was applied again. 
This process continued until the terms indicating a 
variation with angle of attack were no longer signif- 
icant to MSR, In areas of highly nonlinear aerody- 
namics, this bin width was 0.5°; in the more benign 
a ranges, a bin width of 1.0° was adequate. Be- 
cause of the amount of maneuvering in the range of 
a from 24° to 28° during each flight, there were al- 
ways enough data to make such small bins amenable 
to analysis. 

With the bin widths small enough to preclude 
terms indicative of the variation of a parameter with 
angle of attack, nonlinearities in other variables such 
as sideslip or roll rate should become visible (ref. 13). 
An example of such a nonlinearity is demonstrated 
in figure 10 where the effect of the entry of /3 3 
into the equation for the rolling-moment coefficient 
is shown. The variable /3 3 was a consistent entry 
for the two unaugmented flights HD8 and HD 19 
in the bins for 24° < a < 26°. For these bins, 
—4° < 0 < 9°. A similar variation was found in 
data from the RAE (Farnborough) 16.5-Foot Tunnel 
as plotted in figure 10 for the purpose of comparison. 
Although the effect is small for —5° < 0 < 5°, it 
means a reduction in the rolling- moment coefficient 
by almost 70 percent for « 8° and a change in sign 
for the dihedral effect from negative to positive 
at 0 « 5°. Since the flight was asymmetric, no data 
were available for / 3 < —4°. 

Although the wind-tunnel tests indicated a 
(pb/ 2V) 3 dependence for C/ (ref. 6), only two bins 
from the unaugmented flights and no bins from the 
augmented flights required nonlinear terms inpb/2V . 


It is postulated that nonlinear roll rate effects were 
not seen in the flight data because most in-flight 
nondimensional roll rates were less than 0.04, the 
value at which the wind-tunnel nonlinearities became 
apparent. The remainder of the model structure was 
linear in 0,p, r, 8^, and 6 r when these variables were 
significantly excited. In particular, models derived 
from partitioned data from the 1986 series consis- 
tently contained 0 and 6^ in the side-force equa- 
tion, 0,p, and 6^ in the rolling- moment equation, and 
0 and 6fi in the yawing-moment equation. Several 
terms that are normally identifiable in system identi- 
fication applications (particularly 6 r for the side- force 
coefficient and r and 6 r for the yawing-moment coef- 
ficient) were not identifiable from the data analyzed. 

Longitudinal . Both the individual maneuvers 
and the partitioned data were analyzed for model 
structure using MSR. Several of the 1983 flights indi- 
cated pitching-moment-coefficient dependence on roll 
rate p and sideslip through /3 2 . Two of the individual 
maneuvers in the 1986 flights required an a 2 term to 
achieve an adequate model for the vertical-force co- 
efficient. None of the partitioned data gave rise to 
nonlinear a dependence since the bin widths were 
chosen specifically to eliminate such dependence. In 
addition, Cm 6h was not identifiable (did not provide 
a statistically significant improvement in fit) in sev- 
eral of the 1983 individual maneuvers. This lack of 
identifiability is unusual for C mf)h and no reason was 
found for it. 

Derivative Values 

All derivative comparisons were made using the 
adequate model structure as determined by MSR. 
This means that the /? 3 term with is in the 

rolling-moment equation. The effect of adding the 
/? 3 term is to make the linear derivative of 
the nonlinear model slightly more negative when 
compared with the linear model. 

Lateral . The sideslip derivatives are given in 
figure 11 for both the 1983 and 1986 flights. The 
flight-estimated values are compared with published 
wind-tunnel results. The change in side-force coef- 
ficient due to change in sideslip Cy^ is considerably 
less when estimated from flight tests than when es- 
timated from wind-tunnel tests. This result is seen 
consistently in the 1983 series for both canard set- 
tings and in the 1986 series for 6 C = 0°. 

The derivative from the 1983 series does not 
reflect the canard effect seen in the wind tunnel. 
The flight results for 6 C — -10° show approximately 
60 percent of the dihedral effect that was measured in 
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the wind tunnel for 18° < a < 26°. The flight value 
stays negative at a = 30° and the wind-tunnel value 
goes to zero. The flight-derived values for 6 C — 0° 
also generally indicate reduced dihedral effect when 
compared with the wind-tunnel measurements. The 
1986 tests tend to confirm the 1983 results for 8 C = 0° 
except for two values from the unaugmented flight 
HD8. No reason could be found for this discrepancy 
except to note that HD8 was the 1986 flight with 
the largest (3 excitation and in which C^ 3 became 

identifiable. The 1983 series yields C n0 values that 
follow the wind-tunnel trends. However, the 1983 
data for 8 C = -10° indicate C Uf3 < 0 by a = 19°, 
whereas the wind tunnel indicates C n ^ = 0 at a = 
20.5°. The 1983 data for 8 C = 0° indicate C U(i to be 
more negative than that indicated by the wind tunnel 
for 25° < a < 30°. This latter statement is confirmed 
by the 1986 data for 24° < a < 28°. Combining 
directional stability with the dihedral effect C\^ 
through the effect of the inertia ratio Izjlx gives the 
parameter 

J 7 

^V.dyn = Cn 0 ~ J^ C b &mOL 

In highly dynamic maneuvers, the sign of C n/idyn 
indicates stability better than simple C n/r The 
C n ^ dyn for the flight analysis was calculated using 
C n(3 values from the flight data analysis shown in 
figure 11 and approximating Ci Q values from the local 
slopes of the C\ versus (3 nonlinear model curve in 
figure 10. Figure 12 shows that for the 1986 results, 
the /3 3 terms in the rolling- moment coefficient give 
Cng dyn > 0 for 0° < (3 < 5° and a % 25°, but 
its value rapidly decreases to a value of —0.64 for 
/ 3 = ±7°. Reference 3 indicates that C f n ^ dyn > 0 for 
a < 30°. 

Of the three yaw rate derivatives, only the cross 
derivative C{ r was consistently identifiable and still 
showed wide scatter (fig. 13). For the 6 C = — 10° case 
in the 1983 flights, the flight estimates agree well with 
the wind-tunnel values. However, for 6 C = 0° in 1983, 
the C\ flight estimates generally exceed the wind- 
tunnel values for 24° < a < 28°. The 1986 estimates 
(6 C = 0° and —5°) appear to confirm the 1983 flight 
results {6 C = 0°). The sparse flight estimates for Cy r 
and Cn r are presented for completeness. The yaw 
rate and yawing moment were not well-excited in the 
1986 flights because there were no rudder inputs de- 
signed and executed for directional excitation. More- 
over, for the augmented flights (HAD9 and HAD16), 
control inputs were actually used to suppress lateral- 
directional excitation. 


The roll rate derivatives are given in figure 14. 
Both rotary and oscillatory wind-tunnel measure- 
ments are presented for comparison. These tunnel 
measurements are for 6 C = —10° in all figures. The 
discrepancy between sets of values from the oscilla- 
tory and rotary tests has already been noted in ref- 
erence 6. The flight-derived values indicate a loss 
of roll damping for 21° < a < 25° for 8 C = —10°, 
whereas the rotary data indicate such a loss only for 
23° < a < 25°. The oscillatory data indicate no loss 
of damping in this region. The 1983 data for 8 C = 0° 
show damping to be maintained for 22° < a < 28°, 
whereas the 1986 data indicate about one-half the 
level seen in the 1983 data for 6 C = 0°. That dif- 
ference could be due to the smaller excitation of 
roll rate in the 1986 series combined with the non- 
linearities in pb/2V at higher roll rates. The C Up 
values were not identifiable for the 1986 series. 

The control effectiveness of a differential horizon- 
tal tail was usually identifiable for all three coef- 
ficients Cy,C/, and C n . The results are given in 
figure 15. Only Cy 6d from the 1986 series shows 
any consistent difference from the wind-tunnel val- 
ues. The limited number of points plotted for the 
1986 flights stems from the fact that lateral controls 
were used only on the augmented flights (HAD9 and 
HAD 16) and then only in a 5-sec open- loop sequence 
to eliminate excessive bank angle. Therefore, there 
was simply not enough data outside the main area of 
interest in angle of attack. 

Longitudinal . The two primary stability deriva- 
tives Ci a and C ma are given in figure 16. All flights, 
except HAD 16, were flown with the c.g. location at 
12.5-percent c. All wind-tunnel values are referenced 
to that location. Flight HAD 16 was flown with the 
c.g. aft at 22.5-percent c. This led to approximately 
10° less horizontal stabilator for trim in the angle-of- 
attack range from 16° to 40°. A cluster of Ci Q data 
points is shown scattered around the wind-tunnel val- 
ues, and the trends seen in the wind tunnel near 
a = 22° are not apparent from flight. The effect of 
the c.g. position is not apparent in figure 16(b) since 
there are too few data points to verify wind-tunnel 
results. Values for the damping-in-pitch derivative 
Cmy are plotted in figure 17. Pitch rate was well- 
excited in the 1983 series because of cross coupling 
with the large lateral motion. Hence, several individ- 
ual maneuvers and bins provided estimates of C mq 
(fig. 17(b)). However, for the 1986 series, as men- 
tioned earlier, the SAS prevented excessive longitu- 
dinal response and the lateral oscillations were not as 
severe as in the 1983 series. Thus, only four maneu- 
vers out of the entire 1986 series produced an identifi- 
able C mq , and these values are plotted in figure 17(a). 
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The horizontal tail effectiveness C msh is plotted in fig- 
ure 18. This derivative was consistently identifiable 
from the positioned data in both the 1983 and 1986 
series. The average flight value indicates that the ef- 
fectiveness of the horizontal stabilator is degraded in 
flight to about 65 percent of its wind-tunnel value. 

Eigenstructure Analysis of Lateral Modes 

The wind-tunnel-determined roll, Dutch roll, lat- 
eral phugoid, and spiral modes were discussed in ref- 
erence 4 to give a better physical overview of the 
expected HIRM flight characteristics. Figure 24 in 
that report is repeated here as figure 19. Figure 20 
presents damping and frequency as given by the av- 
erage of the flight-derived derivative estimates for 
each angle of attack. Also, flight air density, inertias, 
and mass were used to derive the values in figure 20. 
Significant differences between figures 19 and 20 in- 
clude a change in notation and hp become 

Ari, k 2 , and &3, respectively) and a change in scale for 
the middle graph (&£ to £3). The change in nota- 
tion eliminates reference to roll, Dutch roll, and spi- 
ral modes for the flight-derived values, although the 
eigenvalues from flight-derived derivatives appear to 
separate out these usual lateral modes. A study of 
the eigenvectors (fig. 21) indicates that a different 
set of three modes is present. Figure 21 gives plots 
of eigenvectors corresponding to each flight-derived 
eigenvalue at intervals of 2° in angle of attack. The 
first column indicates a large bank-angle component 
(0) accompanying the heading component ($). The 
eigenvalue that would normally represent roll damp- 
ing is associated with the eigenvector plotted in the 
second column. Even though bank angle (0) mostly 
dominates, a small sideslip component (/?) is present 
and grows steadily until a = 28° where it is then 
equal to </>. Finally, the third column of vectors shows 
that the undamped oscillatory mode is basically a 
roll-dominated (bank angle) mode; it only begins to 
resemble a true Dutch roll at a = 30°. 

Verification 

The final stage of identification, as seen in fig- 
ure 6, is model verification. One phase of verification 
has already been completed by making sure that the 
derivative estimates are physically reasonable and 
comparing them with wind-tunnel test results. The 
other verification technique checks the prediction ca- 
pability of the model by integrating the aircraft equa- 
tions of motion for one of the flight maneuvers. A 
batch simulation was used with the math model de- 
termined by fairing a line through and averaging the 


flight results given in this report. A three-degree- 
of-freedom simulation was used with the longitudi- 
nal response variables taken as measured. Figure 22 
shows the results of that simulation and the angle-of- 
attack time history using a maneuver from HD1. The 
frequencies predicted by the model are good and the 
amplitudes of the oscillations exceed those recorded 
in flight. An additional point of view was gained by 
comparing the predicted and measured time histories 
for lateral force and moment coefficients as shown in 
figure 23. 

Conclusions 

Flight data from two series of flights of the Royal 
Aerospace Establishment (RAE) high-incidence re- 
search model (HIRM) have been analyzed. Model 
structure and the corresponding stability and con- 
trol derivatives were determined. These results were 
compared with available wind-tunnel results, and the 
following conclusions are presented: 

1. Partitioning and modified stepwise regression 
can provide stability and control derivatives from 
data where individual maneuvers were not intended 
or were not suitable for system identification. 

2. Nonlinearities in a variable can be eliminated 
by partitioning, thus allowing for nonlinearities in 
other variables to be visible. 

3. Data from a control-surface-initated response 
at several flight conditions are more useful than data 
from self-excited oscillations because data on and 
beyond the roll rate and sideslip boundaries of the 
self-excited motion can be analyzed, thus giving more 
information on those boundaries. 

4. Estimated parameters from the 1986 flights 
confirmed, in general, the results from the 1983 
flights with a symmetric canard deflection of 0° ex- 
cept for the values of damping- in- roll derivatives. 

5. The vehicle exhibits three unusual roll- 
dominated lateral modes near an angle of attack of 
24°. 

6. The linear model on the lateral variables was 
adequate for each bin except those few where cubic 
sideslip dependence ( ft 3 ) was entered into the rolling- 
moment coefficient yielding a corresponding non- 
linear derivative (C^ 3 ). 

NASA Langley Research Center 
Hampton, VA 23665-5225 
August 21, 1989 
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Appendix 


mg 

qS 


ay = Cy 


Airplane Equations of Motion and 
Postulated Models 

The airplane equations of motion are referred to 
the body axes. They are based on the assumption 
that the airplane is a rigid body. The equations have 
the form 

• p+p v2s r 

u = —qw + rv — a sin u H — Cx 

2m 

^ ^ a ■ r 

v — —ru + pw + g cos V sin 0 + ■ O y 

w - —pv + qu + q cos 8 cos 0 + C% 

2m 

Iy - Iz Ixz , I \ I pV 2 Sb 

p = — j qr + —J — (pq + r) + — — C\ 

lx *X 2I X 

h ~ lx . J XZ / 2 _2\ | pV 2 Sc 

q =^— pr+ —( r -p) + C m 

r=!lLjJr m+ iM-ir-v) + ^ C„ 

0 = p + (q + sin 0 + r cos 0) tan 8 
8 — q cos 0 — r sin 0 


For the stepwise regression method, the equations of 
motion can be formulated as 

mg c 

—£ a X = C X 

qS 


mg r 
—^o-Z = P'Z 
qS 


lx_ 

qSb 


Jy_ 

qSc 

h_ 

qSb 



Iy-Iz 

Ixz 

p- 

Ix 

qr - Ix 

■ 

iz-lx 

Ixz 

Q ~ 

Iy 

pr Iy 

- 

I x - Iy 

Ixz 

r — 

Iz 

pq Iz 


(pq + r) 

(r 2 -p 2 ) 

(p - q r ) 


= Q 
= Cm 
- C n 


For small-amplitude maneuvers, the aerodynamic 
coefficients are postulated as functions of the state 
and input variables and their combinations, and they 
are presented as follows: 

1. The longitudinal coefficients Cx,Cz * and C m 
as functions of a, qcj 2V, <5/,, a 2 , (qc/2V)a, b^a, P 2 , 
a/3 2 , a 3 , a 4 , a 5 , a 6 , a 7 , and a 8 . 

2. The lateral coefficients Cy, C { ■ and C n as 
functions of P, pb/2V, rbj 2V, b d , b T , Pa, (pb/2V)a, 
(rb/2V)a, b d a, b r a, Pa 2 , (pb/2V)a 2 , {rb/2V)a 2 , 
b d a 2 , b r a 2 , P 2 , p 3 , P\ P b , (pb/2V) 2 , (pb/2V) 3 , 
(fia, a, a 2 , and a 3 . 

The variables in both model forms represent the 
increments with respect to their trim values. In 
the equation for the pitching- moment coefficient, the 
term a was not explicitly included to avoid identifi- 
ability problems that occur because of the high cor- 
relation between a and q. Relationships between pa- 
rameters in the expressions for C m (a, P, q, 6^) and 
for C m (a, a, P, q, b^) can be found in reference 13. 
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Table I. Summary of Wind-Tunnel Tests and Figures Containing the Results 


Derivative 


Model size 


Test 


RAE tunnel 


Reference 


Figure 


C L a 


0.36 


Static 


No. 2 11 2' by 8-Foot 


16 


Cm a 

Cm 6h 

Cm„ 


0.36 

.36 

.44 


Static 

Static 

Oscillatory 


No. 2 llj- by 8-Foot 
No. 2 11 5- by 8- Foot 
13- by 9-Foot 


3 

3 

4 


16 

18 

17 


C Y, 


3 


°Ym 

C Y t 


0.36 

.36 

.44 


Static 

Static 

Oscillatory 


No. 2 11 5- by 8-Foot 
No. 2 11 2* by 8-Foot 
13- by 9- Foot 


3 

3 

4 


11 

15 

13 


C b 

C kd 

Cl r 

C b 

Cl 


1.00 

.36 

.44 

.44 

.22 


Static 

Static 

Oscillatory 

Oscillatory 

Rotary 


16.5-Foot 

No. 2 11 y by 8- Foot 


5 

3 

4 

6 
6 


10 , 11 

15 

13 

14 
14 


C n0 

1.00 

Static 

Cn 6d 

.36 

Static 

Cn r 

.44 

Oscillatory 

Cn p 

.44 

Oscillatory 

Cn p 

.22 

Rotary 


16.5-Foot 

No. 2 11 5- by 8- Foot 
No. 2 llj- by 8- Foot 
13- by 9-Foot 
13- by 9- Foot 


5 

3 

4 

6 
6 


11 

15 

13 

14 
14 


13 
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(a) a and V versus time. (b) p and (3 versus time. 

Figure 2. Time histories for selected input and response variables for the beginning of one of the unaugmented 
1983 flights. 
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(c) r and q versus time, (d) ax and a y versus time. 
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Figure 2. Continued. 



120 


o 



oooooooooo 

t- CM CO CM i— r- CM 

■ ■■I >i 


6ap 1 3 2 Bap ‘ P° g 


18 


Figure 2. Concluded. 
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(c) r and q versus time. (d) a x and a y versus time. 




Figure 3. Continued. 




Figure 3. Concluded. 
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Figure 4. Continuation of unaugmented 1983 flight showing spontaneous development of lateral oscillation. 
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(c) r and q versus time. (d) a x and a y versus time. 
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Figure 4. Continued. 
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Time, sec Time, sec 

(a) Differential horizontal tail doublet from (b) Control surface motion that initiated 

1983 test series. individual maneuvers for 1986 test 

series. 




Figure 6. General flow chart for system identification analysis. 
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Figure 7. Angle-of-attack time history and partitioning strategy. 
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Figure 8. Cross plots of roll rate versus sideslip for eight bins after partitioning from one of the 1986 flights 
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Figure 8. Concluded. 
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(a) Measured and computed time histories. 
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(b) Residual time histories. 
Figure 9. Results of compatibility check. 
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(c) Input time histories. 
Figure 9. Concluded. 
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Figure 10. Nonlinear dependence of rolling- moment coefficient on angle of sideslip with 6 C — 0°. 
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Figure 11. Sideslip derivatives. 
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Figure 11. Concluded. 




O Wind-tunnel value 

□ Calculated from flight- 
estimated derivatives 


n p,dyn 


1.2 

.8 

.4 

0 

-.4 

-.8 

- 1.2 




□[ 

H D 

□ c 

n 

u 

u 

_□ 

□ 

1 1 1 

i i i i 


-8 


- 4-2 0 2 

P, deg 


6 8 


Figure 12. Variation of C n/jdyn with /? for a = 25°. 


36 



Wind tunnel 



U) 

i 00 “D 
CM ~ 


O O 


CD 

TJ — 
TJ g 

2 | 

.2 2 

r c 

£ i 


\ ° 
o\ < 

\o 


Figure 13. Yaw rate derivatives. 
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Figure 13. Concluded. 
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Figure 14. Roll rate derivatives. 
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Figure 14. Concluded. 




Figure 15. Differential tail derivatives. 




Figure 15. Concluded. 



Figure 16 . Slopes of lift and pitching-moment coefficients. 
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Figure 17. Pitch rate derivatives. Figure 18. Horizontal tail effectiveness. 
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Figure 21. Eigenvectors using derivatives estimated from flight data. 
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Figure 22. Comparison of measured lateral time histories with those computed from mathematical model of 
HIRM. 
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Figure 23. Comparison of measured lateral aerodynamic coefficients with those computed from mathematical 
model of HIRM. 
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